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ABSTRACT 

A  new  set  of  spherical  shell  theories  is  presented  with  differing  degrees  of  physicality,  varying  from  a 
simple  membrane  model  up  to  a  fluid  loaded  shell  theory  that  includes  translational  motion,  rotary  inertia,  and 
transverse  shearing-stress.  Numerical  results  from  these  theories  are  compared  with  exact  results  from 
elastodynamic  theory.  The  motivation  of  this  study  is  to  overcome  the  limitations  of  both  modal  analysis 
approaches  and  of  the  somewhat  more  general  Extended  Boundapi  Condition  (EBC)  method  due  to  Waterman, 
sometimes  called  the  T-matrix  method.  The  spherical  shell  is  studied  with  an  eye  towards  generalizing  the  results 
obtained  to  spheroids.  The  aim  of  the  present  research  is  to  yield  a  better  understanding  of  the  scattering  event  by 
employing  more  general  and  mote  physical  shell  theories. 

INTRODUCTION 

So-called  "shell  theories"  simplify  the  calculations  of  the  motion  of  thin  elastic  shells  by  making 
assumptions  about  the  scatterer  and  its  movements.  We  use  the  standard  assumptions  of  shell  theory  as 
formulated  by  A.  E.  H.  Love  [1]  and  which  are  as  follows:  First  the  thickness  of  a  shell  is  small  compared  with  the 
smallest  radius  of  curvature  of  the  shell;  second  the  displacement  is  small  in  comparison  with  the  shell  thickness; 
third  the  transverse  normal  stress  acting  on  planes  parallel  to  the  shell  middle  surface  is  negligible;  and  Finally  the 
Fibers  of  the  shell  norjpal  to  the  middle  surface  remain  so  after  deformation  and  are  themselves  not  subject  to 
elongation.  These  assumptions  are  used  in  the  development  of  a  shell  theory  for  an  elastic  spherical  shell  in  the 
spirit  ofTimoshenko-Mindlin[2,3]  plate  theory. 

DERIVATION  OF  THE  EQUATIONS  OF  MOTION 

In  spherical  shells  membrane  stresses  (proportional  to  P)  predominate  over  flexural  stresses  (proportional  to 
P1)  where 


P  = 


1  h 
■y/12  a 


(1) 


We  differ  from  the  standard  derivation  for  the  sphere  [4]  by  retaining  all  terms  of  order  /32in  both  the  kinetic  and 
potential  energy  parts  of  the  Lagrangian  and  by  considering  the  resonance  frequencies  for  the  fluid  loaded  case  to 
be  complex.  We  note  that  this  level  of  approximation  will  allow  us  to  include  the  effects  of  rotary  inertia  and  shear 

distortion  in  our  shell  theory,  as  well  as  damping  by  fluid  loading.  The  parameter  p  itself  is  proportional  to  the 
radius  of  gyration  of  a  differential  element  of  the  shell  and  arises  from  integration  through  the  thickness  of  the 
*hell  in  a  radial  direction.  We  will  use  an  implicit  harmonic  time  variation  of  the  form  exp(-i'ror).  We  begin  our 
derivation  by  considering  a  u,v,w  axis  system  on  the  middle  surface  of  a  spherical  shell  of  radius  a  (measured  to 
raid-shell)  with  thickness  h,  as  shown  in  Fig.  1. 
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I _ __ _ I 

Fig.  1.  —  Spherical  shell  showing  coordinates  used.  The  u,  v,  w  -  coordinate  system  is  somewhere  on  the  mid¬ 
surface. 

Lagraneian  variational  analysis 
The  Lagrangian,  L,  is 


L  =  T-V  +  W, 


(2) 


where  T  is  the  kinetic  energy,  V  is  the  potential  energy,  and  W  is  the  work  due  to  the  pressure  at  the  surface. 
The  kinetic  energy  is  given  by 

1  f2*  pm  pkll  . 


T  =  -p,fo  +  *r?Xo  +  x)2™  edxdddtp. 


where  the  surface  displacements  are  taken  to  be  linear  as  in  Timoshenko-Mindlin  plate  theory: 

„  x..  xdw 


«)] 


and 


w.  =  w. 


The  motion  of  the  spherical  shell  is  axisymmetric  since  the  sound  field  is  torsionless.  Thus  there  is  no  motion] 
the  v-direction.  Substitution  of  Eqs.  (4)  and  (5)  into  Eq.  (3)  yields,  after  integration  over  x  and  6, 

_  r*  A5  A3  ,  A5  A3  <?w  A’  A ,A3  . 

r  =  "’■i“<K«o?+T+to  *  '2W+7)"1e*io7*T2X*>  +(u*ta  )w - 

or,  in  terms  of  p. 


T  =  xpjia1  j*[(\.&P*  +  6/J2  +  l)u2  -Q.6p*  +  6p1)u^H^P*  +  (02  +  l)wJ)sinfttf, 


50 


50 


where  the  first  and  last  terms  in  square  brackets  in  Eq.  (7)  are  associated  with  linear  translational  kinetic  < 
and  the  middle  two  terms  are  associated  with  rotational  kinetic  energies  of  an  element  of  the  shelL 
The  potential  energy  of  the  shell  is 


v  =  \  +  ><x  +  0)2  sin  ******* 


where  the  nonvanishing  components  of  the  strain  arc 


1  (  du  ") 
a{dO  ) 


x  (  du_  d 2 
~"^6 


w 

do2 


(9) 


and 


eM  =  -(cot  du  +  w)  + 
a 


x  J 

—  cot  9 \  u - , 

a  l  de ) 


(10) 


and  where  the  nonzero  stress  components  are 


a„  ~ 


~  ,  2  v  «# 

1  -  V 


Uw  +  vC..l 


(ID 


and 


°h  =  +  Vc-^ 


(12) 


where  E  is  Young's  modulus.  By  substitution  the  potential  energy  becomes 


rK  It  rh!2 

’  E  1  f 

Jo  Jo  J-A/2 

1  -  v2  (x  +  a)2[ 

10-)^  -  i  «■)*<»  «<l  -  — )u  -  i  £l  +  »)■ 

a  d0  a  dO2  a  a  d9 


_  ,  n,,,  x 1  x  dw  x.du  X  d7w 

+2  v(cot  0[  (1  +  -)u - — ]  +  w][(l  +  -)— - —  +  tv] 

a  a  d6  a  d9  a  90 


\x  +  a)2  sin  OdxdOtp , 


(13) 


which  after  integration  is 


nEh  r*,,  ou^2  ,  „  ,  du 

=  - - r  ((w  +  ^rr)  +(”  +  ucotO)2 +2v(w  + —~)(w  +  ucolO) 

1  —  \tl  Jo  IQ 


du 


du „ 


1-  vl 


do 


do 


~  W*' =°''' 9<“ ' -  & * ; 21 ' ■<* ' -  £x  W -  f>s" 1 ™ 


(14) 


Terms  in  the  potential  energy  proportional  to  /? 2  are  due  to  bending  stresses. 

And  finally,  the  work  done  by  the  pressure  of  the  surrounding  fluid  on  the  spherical  shell  is  given  by 


(V  =  Ixa7  f  P'\vs\n9d0, 

Jo 


U5) 


where  pt  is  the  pressure  at  the  surface. 

IkLagranrian  density  and  its  equations  of  motion 

A  Lagrangian  density  must  be  used  instead  of  the  Lagjangian  since  the  integration  along  the  polar  angle  is 
intrinsic  to  the  problem.  The  Lagrangian  density  is 
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(16) 


L  =  ;rp>2I[<l  +  6/3J  +  L8p*)u2  ~(6p2  +  3.6p*)u^  +  tf2  +  \.8pi)(~^)2 

+(1  +  /32)vv2 ]sin  8  -  K-E~  {( w  +  ) 2  +  (w  +  ucot  0)2  +  2 v(w  +  -|^)( w  +  wcot  0) 

1  —  v  08  do 

+2ita2pawsin  8 , 

with  corresponding  differential  equations  of  motion 


and 


0  dL  d  dL  d  dL 
du  d8  du9  dt  du, 


0  -dL  d  dL  d  OL  ^  d1  OL  ^  dr  dL 


dw  dd  dwe  dt  dw,  d8dt  dw9 ,  dQ1  dw9t 
where  subscripts  denote  differentiation  of  the  variable  with  respect  to  the  subscript. 
By  substitution  of  Eqs.  (17)  and  (18)  into  (16)  we  obtain 


0  =  (1  +  /)2)|—  +  cot  0-|^- (v  +  cot2  8)u 
od  dd 


d2w  „  d2w 

-  P 2 - T-B-  cot  8 - ; 

00’  H  d6 


(17) 


(18) 


.  dw  a2 


d2u 


+K1  +  v)  +  P\v+cot2  8)]Z1-^{(L8P*+6P2  + l)~~  (1.8/?*  +  3 P2)^} 

dd  Cp  dt  dddt 


(19) 


and 


- P .  =  P7  0  ■ +  2/J2  cot  00  -  [(1  +  v)(l  +  P2)  +  P2  cot2  0)]J| 

+cot0[(2  -  v  +  cot2  0)/?2  -(1  +  v))u-/3J^^-2/32cot0-^-^ 

aa  ad 

+/?2(1  +  v  +  cot2  0)-^-^ -/?’ cot 0(2-  v+cot2 0)^- 2(1  +  v)w 
ad  ad 

+£![_(1.8/)‘  +  3p2)~  ~  (1.8  /?*  +  3/?2)cot  0^ 


(2C 


These  differential  equations  of  motion  (19)  and  (20)  have  solutions  of  the  form 


«(n>  =  i</.(  i-t?2)"’^, 

»«o 


and 


MV  )  =  '£ws.(r}). 
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where  rj  =  cos 6  and  P.Cl)  are  the  Legendre  polynomials  of  the  first  kind  of  order  n.  When  the  differential 
equations  of  motion  (19)  and  (20)  are  expanded  in  terms  of  Eqv  (21)  and  (22).  «f  obtain  a  vet  of  linear  equations 

in  terms  of  (/„  and  Uy ,  whose  deterroinant  must  vanish.  Wc  shall  consider  two  cases:  w  uii  ami  without  fluid 
loading. 

jhr  vacuum  case 

The  vacuum  case  is  the  simpler  problem  that  occurs  when  the  spherical  shell  is  surrounded  by  a  vacuum 
such  that  there  is  no  damping.  Thus,  the  pressure  at  the  surface  vanishes:  pc  -  0.  The  set  of  linear  equations  the 
expansion  coefficients  must  satisfy  are 

0  =  |i2'(l  +  6/i:  +  l.S/?J)-(!  +  +|ii'CV?  r  i.«/?J)-/f 'v-d  -r  cfill,.  (23) 

and 

0  =  -AJi  v  -  3)/J:  -  l.S/J4  +  1  +  v)(q  +  |f2‘(l  +  2/ri  +  1 .  K/3 ' )  —  2(1  +  i>)-/J?vA.l  IT.  (24) 


where  fi  =  coal  cr,  v  =  v’  +  A.  -  I.  and  A.  =  n(n  +1).  In  order  for  Eqs.  (23)  and  (24)  to  be  satisfied 
simultaneously  with  a  non-trivial  solution  the  determinant  of  the  system  must  vanish: 

0  =  Q4(l  +  6/3-  +  1.8/J4)il  +  2/3’  +  1.8/T  )+Q3((3/3?  +  1.804)A„|<  v- 3)/3?  -  1.8/3"  +  1  +  v| 

-(2(1  +  v)  +  /?jkA.  |(1  +  6/?3  +  1. 8/?4 )  —  (1  +  P2)k(\  +  2pz  +  1.8/?J  )| 

+(1  +/3:)vf2(l  +  v)  +  /J}  vA„  1  -  A.| ( v -  3)/J:  -  1.8/3"  +  1  +  v-)(/J:  +  1  +  v).  (25) 

Since  there  are  no  damping  terms,  the  shell  vibrates  theoretically  forever.  Thus,  the  normalized  frequency 
Q  can  be  taken  to  be  real.  Equation  (25)  is  quadratic  in  Q\  thus  we  expect  two  real  roots  to  (25)  and  thus  two 
modes  for  the  motion  of  the  shell.  They  are  the  symmetric  and  antisymmetric  Lamb  modes. 

The  fluid  loaded  case 

The  fluid  loaded  case  requires  that  we  consider  a  modal  expansion  of  the  surface  pressure  in  terms  of  the 
specific  acoustic  impedance  z,.  In  its  most  general  form  this  is 


p(a,6,<?)  =  XX-’  .)V„P.m(cos6)cosm0, 

<i*0  «= 0 


(26) 


where 


It  Aka) 

r,  =  tp^— - . 

h’Jka) 


The  specific  acoustic  impedance  r,  can  be  split  into  real  and  imaginary  pans: 


where 


(27) 


(28) 


r  =  pc  Re 


<7».(Aa)| 
hpka) )' 


(29) 
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and 


"  o)  \  h'(ka) 


For  the  case  of  axisymmetric  motion  we  are  considering,  the  surface  pressure  is 


or  by  substitution, 


P.(0)  =  -2>X'>.(cos0), 


P.(0)  =  -I(-W,  -  co2W.m.)Pm(cos  Q). 


Use  of  Eq.  (32)  in  our  set  of  differential  equations  of  motion  (19)  and  (20)  yields  the  following  set  of  linear 
equations  for  the  expansion  coefficients  in  the  case  of  a  fluid  loaded  spherical  shell: 

0  =  [£22(1  +  6p2  + 1.8/3*)  -  (1  +  P2)k]U,  +  [£22(3/J2  +  1.8J34)  -  P2  k  -  (1  +  v))W„ 


0  =  -A.[( k- 3)P2  -1.8p4  +  l+  v][/.+[£22(l  +  a  +  2p2  + 1.8/?“)- 2(1  +  v)  +  Qiy - p2 kX JH'„,  (34 


Again  the  determinant  of  Eqs.  (33)  and  (34)  must  vanish.  However,  in  this  instance  the  value  of  Q  must  be  t 
to  be  complex;  the  resonances  have  a  width  that  depends  on  the  damping.  The  result  of  setting  this  detenninai 
zero  is  i 

0  =  f24(l  +  6  P2  + 1.8  P*)(l  +  a  +  2p2  + 1.8  P*)  % 

+n3ty(l  +  6p2  + 1.8 p* )  +n2{(3 P2  + 1.8 p*)X'[( k - 3)P2  -\.8p4  +  l  +  v]  j 

-[2(1  +  v)  +  /?2kA„](1  +  6  P2  + 1. 8  p* )  -  (1  +  P2)k(  1  +  a  +  2P2+ 1.8/?“)}  +f2(-ty(l  +P2)k]  J 

+(1  +  p2)K[ 2(1  +  v)  +  P2kX,]~  A.[(»c - 3)P2  -\.8p4+\  +  v]{p2  + 1  +  v).  J 

Equation  (37)  has  at  least  four  complex  roots.  From  work  with  an  exact  modal  solution  to  the  probtai 
expect  two  roots  to  be  associated  with  the  symmetric  and  antisymmetric  modes  of  the  shell.  We  expect  thy 
two  roots  to  be  associated  with  a  water-borne  pseudo- Stoneley  wave.  j| 

CONCLUSIONS  a 

The  next  step  is  to  plot  the  roots  of  Eqs.  (25)  and  (37)  to  compare  the  resonances  predicted  by  these  3 
with  those  given  by  exact  modal  expansion  solutions.  By  suppressing  a  and  y,  the  model  associated  wjJ 

(37)  reverts  to  the  vacuum  case  model  associated  with  Eq.  (25).  Similarly  suppression  of  factors  of  P  in  Bfl 
will  result  in  a  reversion  to  a  previously  derived  solution  (Junger  and  Feit,  1 986).  We  may  then  rank  tbfl 
different  models  according  to  their  degree  of  physicality  and  compare  their  results  for  various  relatiT® 
thicknesses  against  each  other  and  against  the  exact  results  of  the  modal  expansion  model.  We  may  alsoCjM 
the  limitations  of  each  of  the  models  including  the  exact  solution,  as  well  as  those  of  shell  models  in  genenj® 
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By  setting  the  values  of  a  and  y  in  Eq.  (37)  to  zero,  we  revert  the  shell  theory  model  to  one  without  fluid 

loading.  Similarly,  by  setting  /?  to  zero  as  well,  the  model  reverts  to  a  membrane  model.  These  models,  fluid 
loaded,  vacuo  case,  and  membrane,  are  successively  less  physically  sophisticated  and  give  successively  less  good 
comparison  with  exact  (modal  expansion)  results.  Starting  with  the  least  sophisticated  model,  we  see  in  Fig.  2 
thick  spherical  steel  shell  dilatational  (symmetric)  and  flexural  (antisymmetric)  mode  resonances  calculated  by  the 
membrane  model.  Here  and  in  the  succeeding  figures  thick  means  lila  =0.1;  thin  means  h/a=0.01.  The  shell 

material  is  a  generic  steel  with  density  p,  =  7.7  times  that  of  water,  shear  velocity  vJ  =  3.24  km/s,  and  longitudinal 

velocity  v,  =  5.95  km/s.  The  surrounding  fluid  is  taken  to  be  water  with  density  p  =  1000  kg/m^  and  sound 

velocityc  =  1.4825  km/s.  The  symmetric  mode  shows  a  good  comparison  between  exact  and  shell  theory 
dictions,  ^ul  antisymmetric  shell  theory  results  for  this  approximation  compare  poorly  with  the  exact 
flexural  results.  Note  that  some  symmetric  mode  resonances  were  not  found  by  our  exact  theory  algorithm.  In 
pjg  3  we  see  thin  spherical  steel  shell  dilatational  (symmetric)  and  flexural  (antisymmetric)  mode  resonances 
alculated  by  the  membrane  model.  Again  there  is  good  comparison  between  dilatational  (symmetric)  mode 
resonances  calculated  by  the  two  methods,  except  for  the  first  couple  of  resonances.  Only  a  few  exact  flexural 
resonances  were  picked  up  by  our  algorithm.  And  again  the  shell  theory  flexural  (antisymmetric)  mode 
resonances  show  the  do  not  asymptote  properly  with  increasing  order.  In  Fig.  4  we  have  thick  spherical  steel 
shell  dilatational  (symmetric)  and  flexural  (antisymmetric)  mode  resonances  calculated  by  shell  theory  without 
fluid  loading  (vacuum).  As  in  the  membrane  model  the  shell  theory  and  exact  calculations  compare  well  for  the 
dilatational  (symmetric)  mode  resonances.  In  contrast  with  the  membrane  model,  however,  the  exact  and  shell 
theory  calculations  for  this  model  show  much  better  agreement  for  the  flexural  (antisymmetric)  mode  resonances. 
This  model  does  not  include  fluid  loading,  but  does  include  the  effects  of  shear  distortion  and  rotary  inertia.  The 
vacuum  shell  theory  flexural  mode  resonances  do  not  asymptote  for  large  size  parameter  ka  to  the  exact  results, 
however.  In  Fig.  5  we  see  thin  spherical  steel  shell  dilatational  (symmetric)  and  flexural  (antisymmetric)  mode 
resonances  calculated  by  shell  theory  without  fluid  loading  (vacuum).  As  in  the  membrane  model  the  shell  theory 
and  exact  calculations  compare  well  for  the  dilatational  (symmetric)  mode  resonances  except  for  the  first  couple  of 
resonances.  This  vacuum  model  does  not  have  fluid  loading,  and  has  insufficient  damping  for  the  first  two 
dilatational  (symmetric)  mode  resonances.  Again,  the  flexural  (symmetric)  mode  resonances  show  roughly  the 
correct  behavior,  but  it  is  not  possible  to  tell  what  the  asymptotic  value  of  the  phase  velocity  would  be  for  large 
size  parameter  on  this  scale.  Next  in  Fig.  6  we  have  a  plot  of  thick  spherical  steel  shell  dilatational  (symmetric) 
and  flexural  (antisymmetric)  mode  resonances  calculated  by  shell  theory  with  fluid  loading.  As  in  the  vacuum 
case  as  well  as  for  the  membrane  model,  the  dilatational  (symmetric)  mode  resonances  compare  well  for  exact  and 
shell  theory  methods.  The  flexural  (antisymmetric)  mode  resonances,  as  calculated  by  shell  theory  with  fluid 
loading,  do  not  appear  to  have  the  correct  asymptotic  limit  for  large  size  parameter,  although  they  do  exhibit 
roughly’ the  correct  behavior  for  lower  values  of  ka.  Finally,  in  Fig.  7  we  see  thin  spherical  steel  shell  dilatational 
(symmetric)  and  flexural  (antisymmetric)  mode  resonances  calculated  by  shell  theory  with  fluid  loading.  The 
exact  and  shell  theory  calculations  agree  well  for  the  dilatational  (symmetric)  resonances  and  exhibit  a  marked 
improvement  for  the  First  several  shell  theory  symmetric  mode  resonances.  This  is  due  to  the  inclusion  of  fluid 
loading  in  the  model.  The  flexural  (antisymmetric)  mode  resonances  show  the  appropriate  behavior  on  this  rather 
limited  size  parameter  scale. 
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